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A self-contained theory of the domain wall dynamics in ferromagnets under finite electric current is 
presented. The current is shown to have two effects; one is momentum transfer, which is proportional 
to the charge current and wall resistivity (p w ), and the other is spin transfer, proportional to spin 
current. For thick walls, as in metallic wires, the latter dominates and the threshold current for 
wall motion is determined by the hard-axis magnetic anisotropy, except for the case of very strong 
pinning. For thin walls, as in nanocontacts and magnetic semiconductors, the momentum-transfer 
effect dominates, and the threshold current is proportional to Vo/pw, Vb being the pinning potential. 
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Manipulation of magnetization and magnetic domain 
wall by use of electric current is of special interest 
recently 0- C3- [J S E3] - from the viewpoint of applica- 
tion to spintronics, e.g., novel magnetic devices where 
the information is written electrically, and also as a basic 
physics in that it involves fascinating angular momentum 
dynamics. 

Current-driven motion of a domain wall was studied in 
a series of pioneering works by Berger [jlBEl- In 1984, 
he argued that the electric current exerts a force on the 
domain wall via the exchange coupling 8]. Later in 1992, 
he discussed that a spin-polarized current (spin current) 
exerts a torque on the wall magnetization, and studied 
the wall motion due to a pulsed spin-polarized current || . 
These theoretical works are, however, based on his deep 
physical insight, and seems to lack transparency as a self- 
contained theory. Also, their phenomenological character 
makes the limit of applicability unclear. In view of recent 
precise experiments 0, 0, 0], a general theory starting 
from a microscopic description is now needed. 

In this paper, we reformulate the problem of domain 
wall dynamics in the presence of electric current, and 
explore some new features such as current-induced de- 
pinning of the wall. We start from a microscopic Hamil- 
tonian with an exchange interaction between conduction 
electrons and spins of a domain wall E3 . With a key ob- 
servation that the position X and polarization 4>q of the 
wall are the proper collective coordinates to describe 
its dynamics, it follows straightforwardly that the elec- 
tric current affects the wall motion in two different ways, 
in agreement with Berger's observation. The first is as a 
force on X , or momentum transfer, due to the reflection 
of conduction electrons. This effect is proportional to the 
charge current and wall resistance, and hence negligible 
except for very thin walls. The other is as a spin torque 
(a force on c/jq), arising when an electron passes through 
the wall. Nowadays it is also called as spin transfer |2j 
between electrons and wall magnetization. This effect is 
the dominant one for thick walls where the spin of the 



electron follows the magnetization adiabatically. 

The motion of a domain wall under a steady current is 
studied in two limiting cases. In the adiabatic case, we 
show that even without a pinning force, there is a thresh- 
old spin current, j™, below which the wall does not move. 
This threshold is proportional to K±, the hard- axis mag- 
netic anisotropy. Underlying this is that the angular mo- 
mentum transferred from the electron can be carried by 
both translational motion (X) and polarization (</>o) of 
the wall, and the latter can completely absorb the spin 
transfer if the spin current is small, j s < j". The pinning 
potential Vb for the wall position (X) affects only if it 
is very strong, Vb > K±/a, where a is the damping pa- 
rameter in the Landau-Lifshits-Gilbert equation. In most 
real systems with small a, the threshold would thus be 
determined by K±. Therefore, the critical current for the 
adiabatic wall will be controllable by the sample shape 
and, in particular, by the thickness of the film, and does 
not suffer very much from pinning arising from sample 
irregularities. This would be a great advantage in appli- 
cation. The wall velocity after depinning is found to be 

(*>« vU/jf) 2 -!- ' 

In the case of thin wall, the wall is driven by the mo- 
mentum transfer, which is proportional to the charge cur- 
rent j and wall resistivity p w . The critical current density 
in this case is given by j cr cx V /p w . 

We consider a ferromagnet consisting of localized spins 
S and conduction electrons. The spins are assumed to 
have an easy z-axis and a hard y-axis. In the contin- 
uum approximation, the spin part is described by the 
Lagrangian 011! El 
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hS<j>(cos9-l)- V pin [0] (1) 
| J((V6») 2 + sin 2 0(V0) 2 ) + sin 2 6 (K + K± sin 2 </>)} 



where a is the lattice constant, and we put S(x) — 
S (sin# cos <f), sin 9 sin <f>, cos 0), and J represents the ex- 
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change coupling between localized spins. The longitudi- 
nal (K) and transverse (K±) anisotropy constants incor- 
porate the effect of demagnetizing field. The constants J, 
K and K± are all positive. The term Vp- m represents pin- 
ning due to additional localized anisotropy energy. The 
exchange interaction between localized spins and conduc- 
tion electrons is given by 
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J d 3 xS(x) ■ (c t crc) x 
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where A and c (cr) are the energy splitting and annihila- 
tion (creation) operator of conduction electrons, respec- 
tively, and <r is a Pauli-matrix vector. The electron part 
is given by H c \ = e k c^ k c k with e k = h 2 k 2 /2m. In the 
absence of V p i n and _ffi n t, the spin part has a static do- 
main wall of width A = {J/K) 1 / 2 as a classical solution. 
We consider a wire with width smaller than A, and treat 
the spin configuration as uniform in the yz-plane, per- 
pendicular to the wire direction x. The solution centered 
at x = X is given by 9 = 9 (x — X), tj) — 0, where 
cos9o(x)= tanh(x/A), and sin^o^) = (cosh(x/A))~ . 
To describe the dynamics of the domain wall, it is cru- 
cial to observe that the weighted average of 0, defined 
by o (t) = J (dx/2X)(f>(x,t)sin 2 9 Q (x ~ X(t)) plays the 
role of momentum conjugate to X, and hence must be 
treated as dynamical Neglecting spin- wave excita- 
tions, we obtain the Lagrangian for X(t) and <f>o(t) as 
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X0 O - -K ± NS 2 sin 2 O - V pin (X), (3) 



where V p i n (X) is a pinning potential for X, and N = 
2 AX/ a 3 is the number of spins in the wall. (A is the 
cross-sectional area.) The equations of motion, derived 
from the Lagrangian, L$ — i?i n t, are given by 



(4) 

sin 20o + Tei*, (5) 




J d 3 x\7 x S {x - X) ■ n{x) 



J d 3 xS a (x - X) x n(x). 
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Here So denotes S(x) with 9 — 9q(x — X), = 0o> and 
rip = (cV^c) (/j = x, y, z) is (twice) the spin density 
of conduction electrons. F e \ represents a force acting on 
the wall, or momentum transfer, due to the electron flow, 
while T e i is a spin torque, or spin transfer, which comes 
from the directional mismatch between wall magnetiza- 
tion Sq(x — X) and n(x). We have added a damping 



term (a), which represents a standard damping torque 
(Gilbert damping), T damp = -§S x 5 1]. Note that the 
spin-transfer effect acts as a source to the wall velocity 
via v el = (X/hNS) T el<z . 

To estimate F c i and v e i, we calculate spin po- 
larization n(x) in the presence of a domain wall 
by use of a local gauge transformation in spin 
space |l5|. c{x) = U(x)a(x), where a(x) is the 2- 
component electron operator in the rotated frame, and 
U(x) = m(x) • <7 is an SU(2) matrix with m(x) = 
^sin e °( a: ~- y ) cos 0o, sin g o( x ~ x ) sm Q; cos e o(.^-x) ^ ^ r pj ie 

expectation value in the presence of electric current 
is written in terms of the Keldysh Green function 
in the rotated frame. For instance, n x (x) = [(1 — 
cos 9 ) cos 2 0o — l]n x + (1 — cos#o) cos0o sin0 o ri y + 
sin 6*o cos 0o h z , where n^(x) = —iTr(G^ x (t,t)a IJ ,), 

G x<r,x'a' (*>*') = i ( a l',a'( t ') a ^A t )) bein g the lesser com- 
ponent of the Keldysh Green function |l6(. After a 
straightforward calculation, we obtain 

TT fi2 A 2, ( 2fc + g)x xf x 
■Tel = -KIT, u q Jk<7 — Cfd(e k+q - a - Cka), 
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to the lowest order in the interaction (with wall) 
u q = ~fdxe-^V x 9 (x) = cosh{ : Xq/2) . The distri- 
bution function ] ka specifies the current-carrying non- 
equilibrium state, and P means taking the principal 
value. As is physically expected, F c \ is proportional to 
the reflection probability of the electron, and hence to 
the wall resistivity, as well as the charge current. In 
fact, adopting the linear-response form, f k(T ~ j (e k a) — 
eE ■ vr(df° /de), as obtained from the Boltzmann equa- 
tion (/°: Fermi distribution function, E: electric field, 
v = fik/m, t: transport relaxation time due to a single 
wall), we can write as F c \ = enjR w in one dimension. 
Here n and j are the electron density and current den- 
sity, respectively, and R w = \ 2=- tzt? (u 2 + + u 2 _ ) is the 
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wall resistance |17| . with ( = (kp + — kp^)/(kp + + kp^) 
and u± = Uk F+ ±k F -- More generally, one can prove the 
relation 0] 



F cl = eN c p w j = enR w IA, 



(10) 



using Kubo formula, where p w = R W A/L is the resistivity 
due to a wall 0|, / = jA, and N c = nLA is the total 
electron number. 

Equations l@J and ©, with Q and l(T71Jl constitute a 
main framework of the present paper. We next go on to 
studying them in the two limiting cases; adiabatic wall 
and adrupt wall. 



We first study the adiabatic limit, which is of interest 
for metallic nanowires, where A 3> kZ 1 . In this limit, we 



^-6(q), and by noting (e k+q .- a - e ka ) q=0 = 



take u 2 

2a A 0, we immediately see from Eq.ljSJ) that F e \ = 0, 
whereas 
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remains finite. The spin transfer in this adiabatic limit 
is thus proportional to spin current flowing in the bulk 
(away from the wall), j s = ^ £ fc k x{fk+ - fk-) {V 
I. 1 being the system volume). In reality, the spin cur- 
rent is controlled only by controlling charge current. 
In the linear-response regime, it is proportional to the 
charge current j as j s — rjj, r\ being a material con- 
stant. This parameter can be written as r\ — J2 a ( a + — 
(7 -)/E Q ( <T + + <T -) f° r a w i re or bulk transport, and 
V = E«(*+ - N -)/J2 a ( N + + for a nanocontact 
and a tunnel junction, where cr^ and N± are band (a) 
and spin (±) resolved electrical conductivity and density 
of states at the Fermi energy, respectively, of a homo- 
geneous ferromagnet. For bulk transport in transition 
metals (such as in wire), ry is expected to be small since s 
electrons dominate the conduction. For tunnel junctions, 
in contrast, it may be large (~ 50% [23), since d electron 
contribution will be dominant because of its large density 
of states. 

As seen from Ea. (ll5ll below, the speed of the stream 
motion of the wall is roughly given by v c \ (except in 
the vicinity of the threshold j CI ). For a lattice constant 
a ~ 1.5A and current density j = 1.2 x 10 12 [A/m 2 ] 
0, we have a 3 j/e ~ 250 [m/s]. This speed is expected 
for strongly spin-polarized materials (77 ~ 1) such as half 
metals. In transition metals, where the transport is dom- 
inated by s electrons, r\ would be a few orders of mag- 
nitudes smaller, say ~ 0.01, which may explain the ob- 
served value ~ 3 [m/s] 6]. Current-driven wall velocity 
may thus be useful in determining polarization 77, which 
is of fundamental importance in spintronics. 

Let us study the wall motion in the absence of pinning, 
-Fpin = 0, by solving the equations of motion, (@J and (0) 
in the adiabatic case (F c \ = 0). The solution with the 
initial condition X = <fio = at t = is obtained as 

ncot (^Jf) = x/l - k 2 coth(7t) + 1 {\k\ < 1X12) 
= \/n 2 - 1 cot(wi) + 1 (|k|>1)(13) 
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where k = 2hv e i/(SK±X), 7 = j+^z^= 

7. For \v c i\ < v CT = SKW/2H (i.e., 



\k\ < 1), cot(aX/A) remains finite as t — ► 00, and the 
wall is not driven to a stream motion but just displaced 
by AX — ^- sin -1 k. In this case, the transferred spin 
is absorbed by (f>o and dissipated through K±, as seen 
from Eq. JSJ , and is not used for the translational motion 
of the wall (X); the wall is apparently "pinned" by the 




FIG. 1: Time-averaged wall velocity as a function of spin 
current, j B , in the weak pinning case (Vo < K±/a). 



transverse anisotropy. Thus even without pinning force, 
the current cannot drive the wall if the associated spin 
current is smaller than the critical value 
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Above this threshold, j s > j" (\k\ > 1), this process 
with K± cannot support the transferred spin and the 
wall begins a stream motion. The wall velocity after 
"depinning" is an oscillating function of time around the 
average value (FigQ 
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Js 2 -(j 8 Cr(1) ) 2 , 
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which is similar to the Walker's solution for the field- 
driven case 0. (The bracket (• • •) means time average.) 

We now introduce a pinning potential, V P i n , and study 
the "true" depinning of the wall by the spin-transfer ef- 
fect in the adiabatic limit. Since spin transfer acts as a 
force on fa, the depinning can be better formulated in 
terms of cf>Q. We consider a quadratic pinning potential 
with a range £; 
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where 0(x) is the Heaviside step function. Then the equa- 
tion for <j)Q reads 



^o(i/ + jicos20o) -v (~sin20 o + -y- 



(1 + a 2 

(17) 

where \i = SK±/h and v = 2V X 2 / £ 2 7iS '. This equation 
describes the motion of a classical particle in a tilted 
washboard potential V with (modified) friction. For 
Vd > f CI (= yuA/2), local minima disappear in V and </>o 
is then "depinned" . Then Eq. indicates that </>o starts 
to drift with average velocity (<f>o) = —v e i/(aX), with os- 
cillating components superposed. The time development 
of X(t) is then obtained as X = (v e \/va){\ — e~ vt / a ) for 
\X\ < The depinning of the wall occurs when X(t) 
reaches £. Thus, the critical spin current for depin- 

crf 1) 

ning is given by j s defined above if the pinning is weak 
(Vq ^5 K±/a), while it is given by 



•cr(2) 
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aVoX, 



(18) 
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if the pinning is strong, Vq > K±/a. Since a is usu- 
ally believed to be small Q, we expect that the critical 
current is mostly determined by K±. This seems to be 
consistent with the observations that the critical current 
is larger for a thinner film 0, Q . It would be interesting 
to carry out measurements on a wire with small K± , such 
as that with a round cross section. 

Let us go on to the opposite limit of abrupt wall, 
A — * 0. As seen from Eq.©, the spin-transfer effect 
vanishes. The pinning-depinning transition is thus de- 
termined by the competition between F e \ — eN c p w j and 
F p i n . It occurs when F c \ — NVo/£, giving the critical 
current density 



2V \ 



NVq = 

£eN c p w ena 3 ^R w A 



(19) 



The average wall velocity after depinning is obtained as 
W = &?§M = (ena 3 /2hS)(IR w )/a. This velocity 
vanishes in the limit, A — > 0, due to the divergence of the 
wall mass, M w = h 2 N/(K ± X 2 ). 

For metallic nanocontacts, where £ ~ A ~ a 20] and 
na 3 ~ 1, experiments indicate that the wall resistance 
can be of the order of h/e 2 = 26 kfl [^j. Thus j cr ~ 
(5 x 10 10 x B c [T])[A/m 2 ], where B c = VqX/^S is the 
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depinning field (/is is Bohr magneton). B c r 
(like in Ref.H}) corresponds to j cr - 5 x 10 7 [A/m 

In conclusion, we have developed a theory of domain 
wall dynamics including the effect of electric current. The 
current is shown to have two effects; spin transfer and 
momentum transfer, as pointed out by Berger. For an 
adiabatic (thick) wall, where the spin-transfer effect due 
to spin current is dominant, there is a threshold spin 
current j" ~ (eX/a 3 )max{K±, oVo} below which the 
wall cannot be driven. This threshold is finite even in 
the absence of pinning potential. The wall motion is 
hence not affected by the uncontrollable pinning arising 
from sample roughness if weak-pinning (Vo < K±/a). In 
turn, wall motion would be easily controlled by the sam- 
ple shape, through the demagnetization field and thus 
K±. The wall velocity after depinning is obtained as 
(X) oc y/ {is) 2 — C?s r ) 2 - I 11 contrast, an abrupt (thin) 
wall is driven by the momentum-transfer effect due to 
charge current, i.e., by reflecting electrons. In this case, 
the depinning current is given in terms of wall resistivity 
p w as j cr oc V Q /p w . 

The two limiting cases considered above are both re- 
alistic. Most metallic wires fabricated by lithography 
are in the adiabatic limit, as is obvious from very small 
value of wall resistivity 23]. In contrast, very thin wall 
is expected to bc formed in metallic magnetic nanocon- 
tacts with a large magnetoresistance (called BMR) [ioj . 



A system of recent interest is magnetic semiconductors 
|24| . where the Fermi wavelength is much longer than 
in metallic systems. As suggested by large magnetore- 
sistance observed recently |2l| , magnetic semiconductors 
would be suitable for precise measurement in the thin 
wall limit. 
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